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Abstract 



This paper reviews several applications of photonic side-bands, used by 
Biittiker and Landauer in their theory of traversal time in tunneling [Phys. 
Rev. Lett. 49, 1739 (1982)], in transport and optics of mesoscopic sys- 
tems. Topics include generalizations of the transmission theory of transport 
to time-dependent situations, optics and transport of mesoscopic systems in 
THz electromagnetic fields, and phase- measurements of photon- assisted tun- 
neling through a quantum dot. 



I. INTRODUCTION 

In 1982 Markus Biittiker and Rolf Landauer published a paper on traversal times in 
tunneling [ [[]], that rekindled the interest in an old topic and served as an inspiration for 
much subsequent research. Several articles in this volume address the recent developments 
in this field, which embraces a wide scope ranging from foundations of quantum theory 
to practical questions concerning ultimate speed limits of nanoelectronic components, a 
very characteristic feature of Rolf Landauer's research. The point of the present article, 
however, is not to participate in that particular discussion, but rather to use a technical 
device presented in the 1982-paper to analyze three other physical systems all of which have 
connections to problems Rolf Landauer has been active in. 

The idea of Biittiker and Landauer was to study the sensitivity of tunneling transmission 
coefficient through a potential barrier to a time-periodic perturbation imposed on the barrier. 
As the frequency of the modulation is varied, the inverse of a cross-over frequency was 
identified as a characteristic time-scale for tunneling. The solution to the time-dependent 
Schrodinger equation in the barrier region was written, following Tien and Gordon [ Q , as 

ij)±{x,t,E) = e ±KX e~ iEt/h exp (-^ sin tot 



*An article to commemorate Rolf Landauer's 70th birthday. 
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g±KXg—iEt/h 



(1) 



where V\ is the amplitude of the harmonic time-modulation, uj is its frequency and k = 
{2m[Vo — E]} 1 / 2 /h with Vq the barrier height. The amplitudes of the side-bands at energies 
E ± nhu are given by the Bessel functions J n . Since the time-modulation was introduced 
just as formal device to probe the energy sensitivity of transmission, Biittiker and Landauer 
were interested in the weak perturbation limit, Vi/hu <C 1. In this paper, on the other hand, 
we ask the question: What are the consequences of Eq. (JTJ) in physical situations where the 
external perturbation is not weak? Several examples will be discussed below. In real physical 
systems interactions (impurities, phonons, electron-electron collisions) must be considered, 
and the interaction terms must be included in the Hamiltonian. It is often convenient to 
formulate the analysis in the language of many-body formalism, and instead of working 
directly with Eq. ([!]), we prefer to use the spectral function, which now reads 



A(p, t, t') — exp —i / dtie(p,ti) 
L Jt> 



(2) 



Note that the time-dependent perturbation requires one to move away from the conventional 
energy representation, A(p,E) = 2tt5(E — e(p)). Further, the above expression is a slight 
generalization of Eq.(l) in that by an appropriate choice of the time-dependent single- 
particle energy e(p, t) several different physical systems can be addressed, e.g., with e(p, t) = 
e(p — qA(t)), where A(i) is the vector potential, a uniform time-dependent electric field may 
be included, while setting e(p, t) = p 2 /2m + V\ cos oot we return to situation discussed in 
Ref . [ [l[ , as is readily verified by doing the integral in Eq. (§) , and expanding the result in 
terms of Bessel functions, in full analog with Eq.([l|). We now proceed to the applications. 



II. TIME-DEPENDENT TRANSPORT IN MESOSCOPIC SYSTEMS 

In order that the consequences of Eqs.(^|,0) be visible in an experiment at least the 
following two conditions must be met. First, the system must maintain a certain degree 
of phase-coherence during a significant part of the transport process. Second, the time- 
variation must be different in different parts of the system, and the particles must be able 
to move between these regions. This brings us in the realm of mesoscopic systems. Now, 
the transport in mesoscopic systems can be analyzed extremely successfully in terms of 
conductance formulas, pioneered by Landauer [|!|, and their subsequent generalizations [ |J. 
Generically, one can express the conductance of a mesoscopic system, coupled by ideal leads 
to external reservoirs, as 

9 = jT(e F ) , (3) 

where T = \t\ 2 is the transmission coefficient, and t is the complex transmission amplitude. 
But here we are addressing a problem going beyond the original formulation: we have an 
external time-dependence acting on the system. Further, it would be desirable to include 
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interactions within the mesoscopic regions: these can be very important due to the small 
number of charge carriers which implies less effective screening. 

Recent years have witnessed a flurry of theoretical papers reporting on generalizations of 
the original scattering-approach to interacting and/or time-dependent systems. It would be 
beyond the present purposes to provide a complete list of references, instead we point to two 
chapters [ |6| , representing complementary views, in a forthcoming volume in Handbook 
of Semiconductors. The formulation we adopt here focuses on the tunneling part of the 
currentQ and the current from the left reservoir to the mesoscopic region (we focus on a 
two-terminal geometry) can be expressed as [ |5| 

Mt) = "¥ lL dtl I ^ ImTr { e "" (t " tl)rL ( e ^^') [G < (t,t 1 ) + f L (e)G r (t,t 1 )}} , (4) 

Here the bold-face entities are matrices in the quantum numbers specifying the states in 
the mesoscopic region; T L is the coupling matrix to the left reservoir (the time-dependence 
may be due to external gates which modify the potential barriers between leads and the 
mesoscopic region) and the Green functions G <,r must be calculated in the presence of the 
coupling to the leads. Typically, one would use the Dyson equation to calculate G r while G < 
requires the use of a quantum kinetic equation, e.g. the Keldysh equation. Thus, Eq.(f|) is a 
formal expression for the time-dependent current, nevertheless it appears to form a suitable 
starting point for further calculations, such as those reported in Refs.f |7|, R Pj. 

Many experiments focus on the average current and it is therefore natural to ask whether 
Eq.(|J) could be simplified in this case. Indeed, one finds 

W)) = - T J ^ lh(e) - f R (e)) ImTr | ^ >- ±L (A(e, *)) J , (5) 

where = + exp((e — fi L / R )/kT)} are distribution functions describing the nonin- 

teracting contacts with electro-chemical potential [Ir/l, and A is an object closely related 



to the retarded Green function [ [T(| . This expression is of the Landauer type: it expresses 
the current as an integral over a weighted density of states times the difference of the two 
contact occupation factors. There is an important distinction, however: the quantity in 
curly brackets does not involve, in general, just the transmission coefficient but rather the 
Green function for the fully interacting system, which must be evaluated in the presence of 
interactions (e.g. electron-electron, electron-phonon, and spin-flip). The derivation of Eq.([5]) 
allows arbitrary interactions in the mesoscopic region, however the energy dependence of the 
coupling matrices to right and left contacts must be proportional to each otherf]. 



x The current measured in the contacts also contains contributions from displacement currents. A 
low- frequency theory for these, as well as long-range Coulomb forces is described in Ref.[|6| 

2 This condition is not too restrictive: in most cases the coupling matrices are assumed to be 
constants (in lack of a detailed model) and the conditions for the validity of Eq.([|) are automatically 
satisfied. 
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Let us now apply these results to a simple example. We consider a single, noninteracting 
state with energy eo in the mesoscopic region under the influence of a harmonically varying 
field with amplitude V\. An explicit solution can readily be written down: 

Combining Eqs.(§) and the imaginary part of @ we find that the current can be written as 
(J(t)) = ^ E / de[f L (e) - f R (e)]T(e - khu)J 2 k (£j , (7) 

where T(e) is the elastic transmission coefficient through the mesocopic system. The result- 
ing low-field low-temperature conductance is then 

9~=J £ T^-k^JtQ;) , (8) 

fe= — OO 

which indeed appears as a natural generalization of the standard conductance Eq.@ to the 
time-dependent situation. These expressions bear a very close mathematical resemblance to 
the results obtained by Tien and Gordon and Tucker [ 11], who found 



(J(Vo,t))= £ Jl(^L)l dc ( Vo -khcu/e) , (9) 



i.e., the rectified current of a system biased with V{t) — Vq + Vi cosut is given as a sum of 
dc-currents I^c evaluated at voltages shifted by integer multiples of photon energies. It is 
important to note, however, that Eq. (|?|) was obtained as an explicit calculation for a simple 
time-dependent resonant level keeping the coupling to equilibrium contacts to all orders, 
while Refs. [ 0, |Tl[ only consider the lowest order coupling between the different parts of the 
system, but do not make other restrictive assumptions. It would be interesting to learn more 
about the precise interrelation of these approaches, in particular because Eq.([|) has been 
very succesful in the analysis of recent experiments on semiconductor superlattices which 
are subjected to strong ac-fields originating from free electron lasers [ [T^, [TB| . |ilj ]. 



III. LINEAR OPTICAL ABSORPTION IN THZ-FIELDS 

Our second application concerns the situation where the system under investigation 
is placed in a strong THz-field (FIR-field) and then its optical absorption coefficient is 
measured with a weak probe field in the near infrared part of the spectrum, corresponding 
to near band-edge absorption in GaAs-based semiconductor systems. As we shall see, the 
photonic side-band structure residing in Eq.(^) leads to interesting, and observable effects. 
The dominant features in the measured optical absorption spectrum derive from two basic 
effects: those due to a modified density of states, and those due to excitons, v.i.z. interaction 
effects. We analyze both of these effects; what makes our discussion somewhat nonstandard 
is that we allow arbitrary strength for the THz-field and thus go beyond x^-effects. 
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FIG. 1. The time averaged GDOS for a 2D-system for a range of THz field intensities, 
parametrized by 7 = e 2 E^ Hz /Ahm*L0 3 = (0.2,0.5,0.8,1.1,1.4,1.7,2.0), as a function of scaled 
energy (e e = TkiS). At low intensities one observes a Stark- like blue-shift of the band edge as well 
as finite absorption within the band gap. The blue-shift is given by ej = "/u, and is physically 
interpreted as the average kinetic energy of a classical charged particle with mass m* placed in an 
oscillating electric field of frequency u and strength i?THz- With increasing intensity side bands 
emerge at e = e g + e/ ± 2hu>. (From Ref.[15]) 



Let us first consider the density of states effects. In a recent paper we have shown [ [15 
that in the noninteracting case the time- dependent absorption coefficient axiyji), where ui 
is the frequency of the probe field, can be calculated from 

^)^p(T lW| ) I (10) 
cnn 

where d is the dipole matrix element, n is the refraction coefficient, and the generalized 
density of states [ |16| is defined in terms of Eq.@, 



p(T,u l ) = -Y,A(k,T,u; l ) , (11) 
n k 

where the Fourier-transform of Eq.(§) is taken with respect to t — f, and T = (t + 1')/2. 
Already this simple result contains some interesting physics. In the dc-limit one recovers the 
standard Franz-Keldysh absorption spectrum, i.e., a finite (but exponentially damped) ab- 
sorption in the gap, and characteristic oscillations in the band region. In the time-dependent 



case, shown in Fig. 1, a dynamical Franz-Keldysh effect occurs [[16], [17]]: The overall absorp- 
tion edge experiences a blue-shift, and additional structure appears, both in the gap region 
and in the continuum. 

Excitonic effects require a generalization of the analysis presented above. Rather than 
just focusing on the density of states, one must evaluate the susceptibility function, x r (^ 5 f) = 
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FIG. 2. Linear optical absorption in a quantum well for a range of THz-intensities, 
parametrized by 7 = e 2 E^ z /Ahm*u> s . The photon energy is hto = 2.5 meV, below the ion- 
ization threshold. For small 7 the main feature is red-shifted, while for increasing 7 a blue-shift, 
due to the dynamical Franz-Keldysh effect, occurs. Also, exciton replicas at ±2hu become visible. 
(From Ref.[19]) 



—id(t — t')([P(t), P(t')]), where P(t) is the polarization. In the noninteracting limit one can 
show that lm.Xo(T, u{) oc p(T,u)i), thus establishing a connection to Fig. 1. The Coulomb 
interaction between the electron and hole can be included by considering the Bethe-Salpeter 
equation [ 18|, which we have generalized to include the THz-field nonperturbatively. In the 
photon side-band language the equation reads 

X r n (K + 2nu) = xS,n( k > + 2nu) 

r r/k 

+ ExS,n-n'(k,^ + 2(n + n'V) J j^V{\k - k'|) X ;,(k, w, + 2n'u) , (12) 

where the Fourier representation is defined via 

X r (t, = E / ^x;(^)e^ (t - t,)+m2 - (t+t,) • (13) 

This equation can be solved via standard numerical methods, and the optical absorption 
coefficient is finally obtained from Imx = Sk I m Xn=o(k, Typical numerical results are 
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shown in Fig. 2. The theoretical predictions can be briefly summarized as follows. If the 
THz frequency is smaller than the frequency corresponding to the ionization energy of the 
exciton (i.e. the energy difference corresponding to Is — > 2p transition), the ac Stark effect 
leads to the red-shift seen at low intensity curves. However, the dynamical Franz-Keldysh 
effect of Fig. 1 leads to a blue shift which eventually overcomes the red-shift, and a net blue- 
shift results. Recent measurements performed at the UCSB Free Electron Laser facility, to 
be fully described elsewhere [ ^U] , are in very good agreement with these predictions. On the 
other hand, if the THz frequency is larger than the ionization thresold, the ac Stark shift 
and dynamical Franz-Keldysh effect work in unison, and a blue shift is always observed, 
both theoretically and experimentally. 



IV. PHASE MEASUREMENT OF PHOTON- ASSISTED TUNNELING 

THROUGH A QUANTUM DOT 

The basic conductance formula, Eq.(|3|), involves the absolute square of the transmission 
amplitude. The question is then: Can one measure the phase of the transmission ampli- 
tude? An affirmative answer was given by the recent groundbreaking experiments of Yacoby 
et al. 
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and Schuster et al. [ [22|]. Their experimental protocol runs as follows: A 
magneto-transport measurement is performed on an Aharonov-Bohm ring with a quantum 
dot fabricated in one of its arms. If the quantum dot supports coherent transport, the 
transmission amplitudes through the two arms interfere. A magnetic field induces a relative 
phase change, 27r$/$ , between the two transmission amplitudes, t and £qd, leading to an 
oscillatory conductance g(B) = (e 2 /h)T(B), with 



T{B) = T<°> + 2Re{t*t QD e 2 ™* / * } + 



(14) 



where $ is the flux threading the ring, $o = hc/e is the flux quantum, and where the ellipsis 
represent higher harmonics due to multiple reflections. In the experiments, an oscillatory 
magnetoconductance of this form was clearly observed thus demonstrating coherent trans- 



mission through the dot [ ^T],f^] Furthermore, controlling the potential on the dot with a 
side-gate voltage, allowed measurement of the phase shift of the transmission amplitude. 
The success of these experiments gave rise to a number of other works which concentrated 
on refining the interpretation of the experimental results 
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Yet, the experiments also 
suggest application to other phase-coherent transport processes. One particular example 
which has been of considerable recent interest, both experimentally [ p!2]j27| -|30| and theo- 
retically [ P^-BBIfl, is photon-assisted tunneling. While photon-assisted tunneling (PAT) 
is intrinsically a coherent phenomenon, existing measurements of PAT are insensitive to 
the phase of the transmitted electrons and do not directly demonstrate coherence in the 
presence of the time- dependent field. We have recently proposed [ |38| a measurement of 



photon-assisted tunneling through a quantum dot in the mesoscopic double-slit geometry 
described above. This is, in essence, a combination of the experiments of Kouwenhoven 
et al. [ p8| , |30| where a microwave modulated side-gate voltage gave rise to photon-assisted 
tunneling through a quantum dot, and the interference experiments of [ |21[ and [ |22"f . 

We focus on transport in the neighborhood of a single Coulomb oscillation peak associ- 
ated with a single nondegenerate electronic level of the quantum dot [ |SB| . The effect of the 
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ac side-gate voltage is described entirely through the time-dependent energy of this level 



e(t) = e <y s ) + V x cos cut , (15) 

i.e., precisely of the type considered in above. Now we also emphasize that the static energy 
of the level e depends on the dc side-gate voltage V s . All other levels on the dot can be 
neglected provided the ac amplitude, Vx, and the photon energy, hu, are small compared to 
the level spacing on the dot. 

In the absence of an ac potential, a suitable model for the transmission amplitude tQD(e) 
through the dot is the Breit-Wigner form, 

tQD<£) = t - t „(K) + *r/2 • < 16 > 

where T = Tl + Tr is the full width at half maximum of the resonance on the dot due to 
tunneling to the left and right leads. Eq. (|T6| ) implies a continuous phase accumulation of 
it in the transmission amplitude as the Coulomb blockade peak is traversed. (Note that the 
Breit-Wigner form is exact for a noninteracting system with T independent of energy.) 

In the dynamic case, the simple Breit-Wigner description must be generalized, and the 
object to evaluate is the S-Matrix element [ 10,|37|. Provided interactions in the leads can be 



neglected, the elastic transmission amplitude tq>D( e ) can be written as the energy conserving 
part of the S'-Matrix between the left lead and the right lead 

lim(e',R\S\e,L)=6(e'-e)t Q »(e) . (17) 



The S'-Matrix is simply related to the retarded Green function of the level on the dot, 



including both tunneling to the leads and the ac potential [ tTfl], and we find [ 38 



tqD(e) = -i^r L T R (A(e,t)) . (18) 

In the spirit of the Breit-Wigner transmission amplitude we can use the noninteracting 
(A(e,t)) given by Eq.(^) in further calculations. At finite temperatures one must compute 
^qd = Jd e (—df /de)tQi){e) where f (e) is the Fermi function, and the final result is 

/ r \ 00 i ? r 

tQD = {- T Tn) E Jk{Vjhu>W[- --(^- e (V s ) - kfuo + i-)] , (19) 

k=— oo " 

where ip' is the derivative of the digamma function, and \i is the chemical potential in the 
leads. 

We emphasize that a conventional conductance measurement would yield information 
only about the time average of the square of the transmission amplitude, and the double- 
slit geometry of Ref. [ PT|P^] is necessary in order to probe the phase. Figure |3| shows the 
computed magnitude of £qd (bottom) and its phase (top), as a function of the level energy 
eo(V s ). As compared to the time- independent case (shown as a dotted line), several features 
are noteworthy. The magnitude of £qd shows photonic side-bands, reminiscent of those 



seen in transmission through a microwave modulated quantum dot [ |28fl. However, there is 
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FIG. 3. Temperature dependence of the phase shift A(f) (top panel) and the square of the 
amplitude (bottom) of ^qd- The level-width is T/2 = 0.1, in terms of which the other parameters 
are V\ = 1.0, u = 1.0, and T = (solid line), 0.1 (dashed line), 0.5 (dash-dotted line). For 
comparison, the T = time- independent results are shown as dots. (From Ref.[39]) 



an important difference from the usual case of photon-assisted tunneling. The amplitude 
of the Aharonov-Bohm oscillation is sensitive only to the time average of the transmission 
amplitude £qd- Hence only elastic transmission through the dot contributes, i. e., the 
net number of photons absorbed from the ac field must be zero. The sideband at say 
e = eo(Vs) — ^ corresponds to a process in which an electron first absorbs a photon to 
become resonant at energy e (K), and subsequently reemits the photon to return to its 
original energy. In Ref.[ 38] we have studied the phase as a function of the strength of the 
time-dependent modulation, and find that it is possible to quench the main transmission 
peak, or change the sign of the slope or the phase at resonance by adjusting the ratio V\/%uj 
to coincide with a zero of the Bessel function Jq. 
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